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Abstract 



We investigate the three-state antiferromagnetic Potts model on a simple 
cubic lattice with a cluster flipping Monte Carlo simulation algorithm in the 
temperature region below the transition into disorder at T c \ . We find both the 
well established broken-sublattice-symmetry (BSS) phase at low temperature 
and a new, rotationally symmetric phase at higher temperature, but below 
T c i. The properties of the second phase and the transition temperature to 
the BSS phase are in disagreement with recent cluster- variation and Monte 
Carlo simulation results, but in agreement with simulations by Kolesik and 
Suzuki. 

05.50.+q, 75.10.Hk 



Typeset using REVTgX 



I. INTRODUCTION: 



Due to their unusual properties caused by highly degenerate ground states, antiferro- 
magnetic Potts models have been attracting much attention. Some recent studies focussed 
particularly on the three-state antiferromagnetic Potts model on a simple cubic lattice. 
Field theoretical calculations, renormalization group theory, Monte Carlo simulations, the 
cluster-variation method!, and the coherent-anomaly methodi have been applied to obtain 
information about its critical behavior. Reference 3 gives a detailed history of developments. 

The q-state Potts model is defined by the following Hamiltonian: 

H = Kj25 Si>s ., (1) 

(m> 

where the spins take on the values Sj = 1,2, ...,q, the sum is over nearest neighbor pairs 
on sites on a lattice, and 5 is the Kronecker delta. The lattice can be split up into two 
sublattices a and b, such that all nearest neighbors of any site belong to the other sublattice. 
K = J/ksT is a dimensionless coupling constant. We set J/ks = 1 and use T = 1/K as a 
dimensionless reduced temperature. 

At low temperatures (T 1) the model with q=3 on a simple cubic lattice is known to 
be in the so-called broken-sublattice-symmetry (BSS) phase: one lattice is mostly occupied 
by a single state, while the spin values on the other sublattice are split equally between the 
remaining two states. The two sublattices are therefore not equivalent. Permutation of the 
three Potts states and the two sublattices leads to a six-fold degeneracy at low temperatures. 

At a temperature of T ci = 1.22605(5) as determined by earlier Monte Carlo studiesB the 
model undergoes an order- disorder phase transition, belonging to the XY model universality 
class.! - ! 

Not much attention had been given to the medium temperature region below T cl , until 
Rosengren and Lapinskas (RL)0, based on their recent cluster-variation calculations, claimed 
the existence of three ordered states of different symmetry: the BSS phase below T C 3 = 
0.7715, a 12-fold degenerate phase between T c3 and T c2 = 0.787, and a permutationally 
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symmetric sublattice (PSS) phase from T c2 to T c \. 

Spin configurations in the different phases can best be described through the set of 
concentrations of Potts states on the two sublattices a and b, 

^ = 4E*.*,o (2) 

where i = 1,2, 3,j = a, b, and L is the linear system size. At very low temperature a typical 
BSS state would therefore have c" « 1, c% w « c|, ~ 0, 4 ~ | ~ c 3 . A typical PSS state 
is postulated to have c" = > §, eg = < eg = c| < ~. 

Subsequently, Kundrotas, Lapinskas, and Rosengren (KLR)i performed standard 
Metropolis Monte Carlo simulations, focusing on a sublattice magnetization 
Mj = |(|Pi — P2I + IP2 — P3I + \pI ~ Pi I) j where j again specifies the two sublattices a and 6. 
If the two sublattices have an equivalent set of concentrations as in the PSS state, then 
the difference between the sublattice magnetizations M = M a — should be equal to zero, 
while for a completely ordered BSS state it should be one. For this reason KLR use M as 
an order parameter to distinguish between the BSS and the PSS phase. 

They then use a method developed by Lee and Kosterlitzi to obtain the difference in 
free energy AF between states of M = and values of ±M corresponding to BSS states 
as a function of temperature and system size, and conclude from their data that there is a 
transition from BSS to PSS at a temperature of T = 0.68(1). 

We follow a different approach, using the previously definedi order paramett 
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^1 = ^(4-4), (3a) 

&=\{cl-c\-cl + c h 2 ). (36) 

In a £i-£ 2 -plane perfect BSS states would lie on a circle around the origin at angles of 
= tan -1 (£2/£i) = utt/3, while PSS states are at = (2n + 1)tt/6 with n = 0, ±1, ±2, .... 
The Fourier component of the angular distribution, 

6 = cos(6tan- 1 (6/6)), (4) 
should then be +1 in the BSS phase and —1 in the PSS phase. 



II. RESULTS 



We performed Monte Carlo simulations on simple cubic lattices with linear sizes of L=8, 
16, 32, and 64 and periodic boundary conditions, using efficient modified cluster-flipping 
algorithms!'!. At least 10 6 lattice updates or Monte Carlo steps (MCS) were performed in 
equilibrium for each temperature and system size. 

Fig. |I] shows the distribution of sublattice concentrations, accumulated into a histogram. 
The locations of the peaks in these distributions correspond to the expectation values of 
concentrations of figure 2 in RL. At low temperatures we can see four peaks with relative 
weights 1,2,2,1 as expected for the BSS phase. If the PSS phase were present at higher 
temperatures, we should be able to see three peaks with weights 2,2,2. However, the system 
never develops into a three-peak distribution. Instead, above T « 0.95 it shows only two 
peaks and a broad distribution between the two peaks. 

This broad distribution can be interpreted as result of a uniform angular distribution 
in 0. As we can see in Fig. |2], this is indeed the case: at low temperatures the probability 
distribution of peaks around values corresponding to the BSS state and has minima for 
PSS states. At higher temperatures the distribution surprisingly becomes uniformly flat, 
making no distinction between BSS and PSS or any intermittent states in that phase. 

Finally we show in Fig. [5] a plot of the thermodynamic average (06) as a function of 
temperature and system size. For all system sizes (06) is positive at low temperatures 
and decreases with increasing temperature. It becomes zero clearly below the transition 
to the disordered phase, which demonstrates a symmetry-breaking transition from the BSS 
state into a phase that is rotationally symmetric in the ^1-^2-plane. The strong finite-size 
dependence of (06) is very unusual. This makes an accurate estimate of the transition 
temperature difficult. However, we find (06) to be effectively zero at T = 1.0 for all system 
sizes investigated. We therefore estimate the transition temperature to lie close to 1.0 in 
disagreement with KLR. Our results do agree with those obtained independently by Kolesik 
and Suzuki!. 
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III. DISCUSSION: 



Our results indicate the following: there is a transition from the BSS state into an 
intermediate state at a temperature T rot lower than the critical temperature T cl . In the 
intermediate region, the distribution of the vector order parameter (£1,^2) is fully rotational 
symmetric, very much like the planar rotation or XY models. There is no region where 
the order parameter corresponds to the PSS phase. This actually agrees with the results 
shown in figure 2 of KLR, where |M|, which should be zero for a PSS state, is significantly 
different from zero for all temperatures below T cl . Fluctuations that would cause (\M\) to 
differ slightly from zero would also be expected to decrease with increasing system sizes. 
This is not the case in figure 2 of KLR. 

KLR's histograms of M also show that between T = 1.0 and T cl there is always a 
continuous range of values of M almost equally accessible to the system, extending all the 
way from BSS to PSS states. Therefore we do not think that their order parameter M is 
suited to distinguish between the two ordered phases and locate the transition temperature. 

Qualitatively our results could be explained in the following way: If one looks at the PSS 
state as a configuration with, i.e., state 1 dominating sublattice a and state 2 dominating 
sublattice b, one can start to "rain down" randomly state 3 spins onto the whole lattice. 
This allows for a higher entropy than the BSS state (where one sublattice is basically "frozen 
out"). The randomness is restricted however, since two state 3 spins do not want to be next 
to each other. The increase in entropy comes therefore at a cost in energy. Our numerical 
results indicate, that this trade-off never favors the PSS over the BSS states, while RL 
reached a different conclusion at this point. 

IV. CONCLUSIONS 

In summary we have found evidence for a rotationally symmetric ordered phase in the 
three-state AF Potts model on the simple cubic lattice in an intermediate temperature range 



5 



between the low temperature BSS phase and the disordered phase above T c \. We estimate 
the transition between the ordered phases to take place around T rot m 1.0, significantly 
higher than previously reported. The exact location and the nature of this transition need 
more investigation. 

We thank E. Domany for stimulating discussions. 
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FIGURES 

FIG. 1. Summary histogram of the sublattice concentrations for linear system size L = 64 
after 2.5 x 10 6 MCS per temperature. The curves for different temperatures have been shifted for 
clarity. 

FIG. 2. Histogram of the angle 4> of the vector order parameter from eqns. (3) for L = 16 
after 1.4 x 10 7 MCS per temperature. The curves for different temperatures have been shifted for 
clarity. 

FIG. 3. Average of the Fourier component as a function of reduced temperature. Filled 
circles are for L = 64, white squares for L = 32, open circles for L = 16, and triangles for L = 8. 
Data points are obtained from averaging over 2.5 x 10 6 ,10 6 ,1.4 x 10 7 , and 3.5 x 10 7 MCS per 
temperature respectively. 
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